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We present preliminary results on the proton spin structure function at zero 
momentum, in the quenched approximation. Our calculation makes use of a nonper- 
turbative means of determining the multiplicative renormalization of the topological 
charge density. 
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I. INTRODUCTION 



Many years after its experimental observation, the proton spin crisis remains a puzzle 
0-0 for QCD to address. The problem may be stated as follows: Consider the singlet axial 
current 

N f 

il = Ys^n^f (!) 

Its on-shell nucleon matrix element has the form 

(P, e\ Jl W, e ) = u(p, e) [Gi(fc 2 ) 7 „76 - G 2 (k 2 )k, l5 ] u(f, e') (2) 

where e, e' label the helicity states and k is the momentum transfer. In a naive wave function 
picture Gi(0) can be interpreted as the fraction of the nucleon spin carried by the quarks. 
Experimental determinations lead to an unexpectedly small value of Gi(0), calling for a 
theoretical explanation in the context of QCD. 

In the chiral limit, j\ is n °t conserved; its anomalous divergence is proportional to the 
topological charge density, 

= -^r^^pM ( 3 ) 

One has 

(p,e\d^l\p',e') = -(p,e\2N f q\p',e') = -2MA{k 2 ) u(p, e)ij 5 u(p', e') (4) 

where A(k 2 ) = G^k 2 ) + G 2 (k 2 )k 2 /M and, therefore, A(0) = d(0) . Finding G x (0) thus 
reduces to a calculation of the proton matrix elements of q(x). 

For a determination of this nonperturbative matrix element, we turn to the lattice. To 
proceed, we need a (lattice-) regularized version of q(x), for example 

1 ±4 

t^paT? [n^n pCT ] (5) 



2 4 x 32tt 2 

pLVp(J=±l 

where is the product of link variables around a plaquette. As with any regularized 
operator, the correct approach to the continuum requires a proper renormalization. A 
consideration of the quantum numbers of q L (x) leads to the following simple form |J: 

q L (x) a 4 Z q (P)q(x) + 0(a 6 ) (6) 

A precise determination of Z q is a prerequisite for the calculation of the hadronic matrix 
element. 

Using the numerical (heating) method described in Refs. ||||, we evaluate Z q (0) without 
any recourse to perturbation theory. Subsequently, an estimate of Gi(0) is obtained by 
measuring the matrix element (p,e\ q L \p',e') on the lattice, using Wilson fermions in the 
quenched approximation. 

We should also mention two earlier investigations of this issue: Ref. @ estimates Gi(0) in 
a quenched calculation; the result quoted there must be corrected by renormalization effects 
of the lattice topological charge density. In Ref. |L0| a simulation of unquenched staggered 
fermions is performed, using the "geometrical" topological charge density. 
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II. DETERMINATION OF Z, 



We determine Z q nonperturbatively by a method described in Refs. |7|||. The method 
relies on the fact that Z q (/3) is produced by short ranged quantum fluctuations, and there- 
fore should not suffer from critical slowing down, unlike "global" quantities, such as the 
topological charge. We start from a configuration Cq which is an approximate minimum 
of the lattice action and carries a definite topological charge Q L0 (typically Q L0 ~ 1). We 
construct ensembles C n made out of independent configurations obtained by heating (i.e. 
updating by a local algorithm) the starting configuration C , at a given /3, for the same 
number n of updating steps, and average the topological charge over C n at fixed n. Plotting 
Q L = 1 L ( x ) averaged over C n as a function of n, we should see first a decrease of the 
signal, originated by the onset of Z q (f3) during thermalization of the short-ranged modes, 
followed by a plateau. The average of Q L over plateau configurations should be equal to 

On a 16 4 lattice we constructed an approximate instanton configuration whose charge 
is Q L ~ 0.95. Heating was performed using a 10-hit Metropolis algorithm (tuned to 50% 
acceptance) at (3 = 6.0, and collecting about 4500 heating trajectories. In Fig. [I] we plot 
Q L (C n )/Q L0 . We see clearly a plateau starting from n ~ 15. The value of Q L /Q L0 at the 
plateau gives an estimate of Z q (f3): 

Zq (f3 = 6.0) = 0.18(1) (7) 

In order to check the stability of the background topological structure of the initial config- 
uration, after the heating procedure we cool the configurations, by locally minimizing the 
action; we find Q L ~ Q L after a few cooling steps. 



III. SPIN CONTENT OF THE PROTON 



In order to calculate Gi(0) we performed Monte Carlo simulations in the quenched 
approximation using Wilson fermions. In terms of the proton annihilation operator a a = 
e abc [u b ^fi'jsd c }u'^ and its Fourier transform a(p,t), we should have 

(S a (0, t)ap(0, 0)) ~ Z a e~ Mt t » 1 (8) 

and 

(M0,t)Uk,t q )M0,0)P, a ) * i(-) p V Z a Z q ^ {k) A{ ~p \k\ e~ Mt e -(*W-*ft t,t q » 1 

(9) 

where E{k) = \/M 2 + k 2 ; P a p and P are the helicity projection operator and its eigenvalue. 
Given that, with the Wilson fermion discretization, the state of opposite parity propagates 
backward in time, we must keep t < T/2, where T is the temporal size of the lattice, in order 
to isolate the proton state. Then Gi(0) = A(0) can be obtained extrapolating to k = the 
following relationship 



3 



^rf 1 = -i(-) P ^^ e ^)-^-^(a Q (0,t)gJ^ga,(0 ) 0)P, Q ) (10) 

when t,t g ^> 1. In practice, since in a finite lattice we cannot reach k — 0, the smallest 
available nonzero momentum = (2n/L, 0, 0) can be used to estimate A(0). 

Our measurements were performed on a sample of about 100 configurations generated 
on a 16 3 x 32 lattice at (3 — 6.0. We computed the Wilson fermion propagator using the 
smearing technique described in Ref. JIT] : we fixed the Coulomb gauge and used a delta 
function localized on a 6 3 cube in the Dirac equation. We considered the following values of 
the hopping parameter: k = 0.153, 0.154, 0.155. 

We found the quantity ( ]T0| ) to be zero within errors, so we can only conclude giving a 
bound on G\(0) (within one standard deviation) 



Gi(0) < 0.025-/ ~ 0.4 (11) 



While the above bound is not stringent enough for a comparison with experimental data, 
the situation is actually quite promising: A higher statistics calculation is called for, and we 
are currently undertaking it on an APE/Quadrics dedicated machine. Also, an improved 
lattice topological charge density operator (less noisy or with Z q closer to one) would be 
welcome. Finally, an obvious further question to address is the effect of unquenched fermions. 
We will be reporting on these issues in future publications. 
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FIGURES 

FIG. 1. Q L /Q L0 versus the number of updatings when heating an instanton configuration at 
(3 = 6. A dashed line indicates the value of Z q estimated by averaging data on the plateau. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9402019vl 



